December 2, 2008 

WEAK AND STRONG TAYLOR METHODS FOR NUMERICAL 
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Abstract. We apply results of Malliavin-Thalmaier-Watanabe for strong and 
weak Taylor expansions of solutions of perturbed stochastic differential equa- 
tions (SDEs). In particular, we work out weight expressions for the Taylor 
coefficients of the expansion. The results are applied to LIBOR market models 
in order to deal with the typical stochastic drift and with stochastic volatility. 
In contrast to other accurate methods like numerical schemes for the full SDE, 
we obtain easily tractable expressions for accurate pricing. In particular, we 
present an easily tractable alternative to "freezing the drift" in LIBOR market 
models, which has an accuracy similar to the full numerical scheme. Numerical 
examples underline the results. 



Let (fi, F, P) be a probability space carrying an ./V-dimensional Brownian motion 
(Wt) t>0 with a d x d correlation matrix. We consider smooth curves F c : M — > 
L 2 (ri;~R N ) of random variables, where e <E R is a parameter. We apply Taylor 
theorems to obtain strong approximations of the curve F e at e = and we apply 
partial integration on Wiener space to obtain weak approximations of the law of F e 
for small values of e. 

We choose the notion Taylor expansion instead of asymptotic expansion in order 
to point out that the strong method is indeed a classical Taylor expansion with usual 
conditions for convergence. The weak method represents a truncated converging 
power series in the parameter e if - for instance - the payoff / : R w — > R stems 
from a real analytic function and some distributional properties are satisfied. 

2. Weak and strong Taylor methods - Structure Theorems 

We introduce in this section two concepts of approximation. Consider a curve 
e i-> F € , where e £ R and F e £ L 2 (Q; R N ). 

Definition 1. A strong Taylor approximation of order n > is a (truncated) 
power series 
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1. Introduction and Setting 



(2.1) 




such that 



(2.2) 



E(\F e -T?(F e )\)=o(e n ) 
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holds true as e — ► 0. 

Remark 1. In our setting a strong Taylor approximation of any order n > of 
the curve F e can always be obtained, see for instance [KM97J. 

Let / : R N — > R be a Lipschitz function with Lipschitz constant K, then we 
obtain 

(2.3) E(|/(F e ) - f{T?(F e )) |) < KE(\\F e - T™(F e )\\) = Ko{e n ). 

Equation (|2.3|) does not hold anymore if / is not globally Lipschitz continuous. 
In particular, we observe the dependence of the right hand side on the Lipschitz 
constant K. Hence, truncating an a-priori known Taylor expansion leads to an 
error term, which contains the Lipschitz constant and is therefore not useful for 
non-Lipschitz claims. The weak method navigates around this feature by partial 
integration. 

Definition 2. A weak Taylor approximation of order n > is a power series 
for each bounded, measurable f : M. N — > M., 

n j 

W?(/,F e ) :=^lE(/(Fo)^), 

i=0 

where 7Tj € L 1 (f2) denote real valued, integrable random variables, such that 

|E(/(F e ))-W?(/,F e )l=°(0- 

Remark 2. The weights iti for i > 1 are called Malliavin weights. 

Remark 3. If the law of F e is real analytic at e = in the weak sense, i.e. if there 
exist (signed) measures \ii such that for all bounded, measurable f : R N — > M the 
following series converges and the equality 

E(f(F e ))=J2^ f f(x)iH(dx), 

i>0 l ' JrN 

holds true, precisely then we do have a converging weak Taylor expansion. We aim 
for constructing stochastic representations of the following type, for i > 0: 

f f(x)f H (dx)=-E(f(F )7r i ). 

For the definition of the weak Taylor approximation to make sense, existence of 
the Malliavin weights has to hold. The following theorem can be found in a slightly 
different version in |MT06j and goes back to S. Watanabe. For the definition and 
notion of T> oc (R N ) see |Mal97] or |Nua06j . 

Theorem 1. Let F e : K -> V°°(R N ) be smooth and assume that the Malliavin co- 
variance matrix j(F e ) is invertible with p- integrable inverse for every p > 1 around 
e = (i.e. on an open interval containing e — 0). Then there is a weak Taylor 
approximation of any order n > and there are explicit formulas for the weights 
7Ti . If we only know that the Malliavin covariance matrix 7(^0) is invertible with 
p-integrable inverse, then we can also calculate the Malliavin weights, since they 
depend only on 7(^0) ■ 
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Proof. Fix n > and take a smooth test function / : M. N — > E and assume that 
7~ 1 (i 7 ' £ ) exists as a smooth curve in T>°° on a open e-interval containing e = 0. By 
standard arguments we can prove the following formula 

|E(/(F e )) =K^(F e )5(a»(D.F e ) T 7-\FJ±F e )\ 

More precisely, by the integration by parts |Nua06| Definition 1.3.1-(1.42)], the 
chain rule [Nua06, Proposition 1.2.3] and the definition of the Malliavin covariance 
matrix, |Nua06| page 92], we obtain from the right hand side the desired left hand 
side. 

Notice that the e-dependence of the Skorohod integral is smooth due to basic 
properties of 1?°°. Hence, we can calculate higher derivatives of the left hand side 
by iterating the above procedure and differentiating the Skorohod integral. Wc 
denote 

(2.4) tt! := S(a » (D s F e ) T ^~ 1 (F e )^F e 



de 

We write then, pars pro toto, the formula for the second derivative 
— E(/(F e )) =Eif(F e )S(s ^ Trxp.FejV 1 ^)^ 

dF T dF 




E 



(f(F e )8(s^(D s F £ f 7 -\F e ) d ^ 7 -\F^ 



+ E\f(F e )5(s^ (D.Fef-f-^Fe) 




This formula makes perfect sense at e = and - by induction - we see that we can 
perform this step for any derivative. The general, recursive result is the following: 

dF 

a s := (D s F e ) T -f^ 1 (F € )—^ for < s < T, 
de 

, d 

7T n := d(s i ^ a s 7r n _i) + — 7r„_i, 
ae 

7T := 1. 

Here we understand the weights ir n as e-dependent, whereas in the final formulas 
we put e — 0. This proves the result for smooth test functions / and under the 
assumption that the Malliavin covariance matrix is invertible around e — 0. If we 
approximate a bounded, measurable function / by smooth test functions we obtain 
the desired assertion by standard arguments, since the weights are integrable. □ 

Remark 4. By Taylor's theorem and the Fad- di- Bruno- formula we obtain 

\ct\<n 

where p a is a well-defined polynomial in derivatives of the curve e m f e , for a 
multi-index a. Since T>°° is an algebra, see iM al97j . the above expression lies in 
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L p (fl) for each p > 0. The previous result provides a representation of the partial 
integration result for 

E (£ fM{F e )p a )=E{f(F e )v n ). 

\a\<n 

The structure of the weights is seen from above. The result can be considered as 
a dual version of the Faa-di-Bruno-formula. However, the structure of this dual 
formula is much simpler. 

We provide an example to demonstrate the strong and weak method of approx- 
imation. The method works in order to replace time-consuming iteration schemes, 
like the Euler-scheme, by simulations of "simple" Ito integrals. 

Example 1. We deal with a generic, real-valued random variable over a one- 
dimensional Gaussian space, see [Nua06] . i.e. 

oo i 

where the F l lie in the (i+ l) st Wiener chaos (one can think of a Hermite 

expansion for instance) and the sum is understood in the L 2 -sense. From the strong 
expansion we obtain immediately - for a given Lipschitz function f : R — * R - that 

|E(/(F e )) -E(/(F° + e.F 1 ))| <Ko(e), 

as e — > 0, where K denotes the Lipschitz constant of f. This simple approximation 
can be sometimes quite useful. 

We assume now that F° = J °° h(s)dW s has non-vanishing variance in order to 
calculate the weights, which do depend only on 7(Fb). The strong Taylor approxi- 
mation is given by definition, the weak Taylor expansion can be constructed by the 
previous recursive formulas and the specifications 

D S F° = h(s), 




Jo°° h (s) 2 ds' 

In order to obtain a first-order approximation for bounded, measurable random 
variables we therefore have to calculate 

E{f(F°))+eE{f(F°)m), 

where 

This amounts to an integration of f times a polynomial with respect to a Gaussian 
density, since: 

/•oo /*oo 

E(/(F°K) =E(f(F°)F 1 / a s dW s )- E(f(F°)D s F 1 a s )ds. 

Jo Jo 

Notice that the strong approximation does not yield such a result for bounded, 
measurable random variables. Notice also that in the given case the approximation 
can be calculated in a deterministic way, since we deal with Gaussian integrations. 
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The second-order weak Taylor approximation is given by 

E(/(F )) + eE(/(^ ) 7 r 1 ) + e 2 E(/(FV 2 ), 

where 

7T2 = <5(s i ► iria s F ) + <5(s i ► a s F ) . 

3. Applications from Financial Mathematics 

For applications we want to deal with strong and weak Taylor approximations 
of a given curve of random variables. We are particulary interested in cases, where 
the first derivative ^P-| e= o is of simple form or - even more important - where the 
Malliavin covariance matrix 7(Fb) is of simple form. In these cases it is easy to 
obtain first or second order approximations of the respective quantities in the weak 
or strong sense. 

In what follows, first we will present one of the most applied interest rate models, 
namely the LIBOR market model (LMM). Then, we will introduce the commonly 
used technique of freezing the drift. We will show how to embed the" freezing the 
drift" technique into our framework of Taylor approximations. We understand 
freezing the drift as a strong Taylor approximation of order zero in the drift term 
of the LIBOR SDE. Our goal is to put this technique into a method, where we 
can in particular improve the order of approximation. We will finally extend the 
assumption of log normality and develop a stochastic volatility LMM, where we will 
show how to obtain tractable option prices via our weak Taylor approximations. 

3.1. The LIBOR Market Model. We apply our concepts to the LMM, initially 
constructed by [BGM97] . [MSS97] and |Jam97| . Let T denote a strictly positive 
fixed time horizon and (£1,J-t,V, (^t)o<t<T) be a complete probability space, sup- 
porting an TV-dimensional Brownian motion Wt = (W*, ...VF t /v )o<t<T- The factors 
are correlated with dW l t dW( = pijdt. Let = T < T x < T 2 < . . . < T N < 
Tjv+i =: T be a discrete tenor structure and a := T^+i — Ti the accrual factor for 
the time period [T^T^+i], i = 0, ...,7V. Let P(t,Ti) denote the value at time t 
of a zero coupon bond with maturity T.- L £ [0, T]. The measure P is the terminal 
forward measure, which corresponds to taking the final bond P(t, T) as numeraire. 
The forward LIBOR rate L\ := L t (Ti,T i+1 ) at time t < T t for the period [T i7 T l+1 ] 
is given by: 

L\ = L t (Ti,T l+1 ) = — ( ' - - 1 
a \P{t,T i+ i) 

We assume that for any maturity Ti there exists a bounded, continuous, determin- 
istic function a l (t) : [0,Tj] — > R, which represents the volatility of the LIBOR L l t , 
i = l,...,N. The log normal LIBOR market model can be expressed under the 
measure P as: 

(3.i) dL\ = a i {t)L\(- J2 'f'ffj Pii) dt + ^mtdwi i = 1, N. 

j=i+l 1 + a ^t 

3.2. Freezing the Drift. The dynamics of forward LIBORs for i = 1, ...,7V — 1 
depend on the stochastic drift term — ^ , i < j < TV, which is determined by LI- 
BOR rates with longer maturities. This random drift prohibits analytic tractability 
when pricing products that depend on more that one LIBOR rate, since there is 
no unifying measure under which all LIBOR rates are simultaneously log normal. 
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In addition, it encumbers the numerical implementation of the model. Common 
practice is to approximate this term by its starting value i ° L ° j or as it is widely 
referred to as freezing the drift, i.e. 

aL\ aL 3 



1 + aL\ 1 + aL 3 Q 

It was first implemented in the original paper [BGM97 for the pricing of swaptions 
based on the LMM. [BWOOj and jSch02j argue that freezing the drift is justified due 
to the fact that this term has small variance. However, by freezing the drift there 
is a difference in option prices with the real and the frozen drift. It has not been 
examined how big the error is or for which assets it works well or not. Our aim is 
to investigate such a phenomenon and improve the performance by providing with 
correction terms of order one. 

3.3. Correcting the Frozen Drift. The purpose of this section is to embed the 
well-known and often applied technique of freezing the drift into the strong and 
weak Taylor approximations, in order to develop a method to improve the order of 
accuracy. Specifically for the strong Taylor approximation, the method works well, 
since we always deal with a globally Lipschitz drift term x <— > with small 

Lipschitz constant a. 

Remark 5. As it will be clear later, the strong Taylor correction method can be 
accommodated with any extension of the l og norm al LMM, for example with the 



Levy LIBOR model by Eberlein and Ozkan EO05) . 



3.3.1. Strong Taylor Approximation. We first state a useful lemma, asserting that 
we can indeed freeze the drift under special model formulation and choice parame- 
ters. 

Lemma 1. Let t\ 6 R and consider for i = l,...,N the following stochastic 
differential equation: 



(3.2) dX^' ei) =ei a'(*)X t (,,ei) ( 



j=i+i 1 + 0iX t ' 1 



defined on the complete probability space (f2, J-r,P, (Tt)o<t<T) where Wt is an N- 
dimensional Brownian motion under the measure P with dW\dWl — Pijdt. Then 
the first-order strong Taylor approximation for X^ 1 '^ is given by: 

(3.3) Ti(X^) = X^+ ei A 



£1=0 ' 



Proof. By ([1} we obtain for n = 1: 

Tl^X^) ~ = + e t Yi + ( ei ), 

since xj: 1 ' ^ — Xq 1 ' ^ and where Y t l :— ^-| ei =o-^t is the first-order correction 
term. By differentiating (|3.2p with respect to t\, we calculate: 

1 J'=»+1 + 
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and derive as the solution to the above linear SDE: 



Jo v jC7ti i + air ' J ° 



with r i = 0. □ 

Remark 6. We parametrise the LIBOR market model in terms of the parameter 
ei as follows: 

and assume at t = i/iai Lq^ 1 ^ = -XTg 1 ' 61 '' /or aZ/ ei and aH i = 1, N. If t\ = 1, 
what we obtain is the standard LIBOR market model formulation and in particular 
Lj* = x[ l,1 \ For e\ = 0, X^ 1 ' ^ equals its starting value and thus the drift term 
in the following SDE is no longer stochastic 



N vti-o) in\ 

(3.5) dl.™ = a\t)L^ ( - £ J(m P- dt + dW i 



3=l +i 1 + aX. 



The next proposition provides a way for a pathwise approximation of , by 

means of adjusting its SDE. This is achieved by adding T™ (X^' ei ^) in the frozen 
drift part. 

Proposition 1. Assume the setup of Lemma\l\ and assume further at t — that 
£V»>«i) _ Xq'^ 1 ^ for all e\ and all i — 1,...,N. Then the stochastic differential 
equation for l[ 1 ' (Li " > with the unfrozen drift: 

(3.6) dLf^ =a\t)L ( i^(- Y aAt ° [t > PlJ dt + dWl 

can be strongly approximated as e\ j, by 

Remark 7. For n = 0, we derive the "freezing the drift" case. For n = 1, we 
already obtain an improvement. 

Proof. First step is to interchange with (xj: J ' )+ in (|3.6|) to obtain: 

V ^ l + a(*<*«>)+ J ^ 
This yields no change for the dynamics of L^ :ei \ since X^ ,ei ^ = (X^' ei ')+. 
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By Taylor's expansion, wc know that as ei J, 0, L[ l ' €l ^ — ► L^ 1 ^ P-a.s. The 
estimate for the error term is given by 

logi^-logL^ 1 ^ 

= / cr (s)(- > — — r^l /»i7+ > , v ' pi, )ds < 

< / t a|Xp^)-(T^ i (Xi^))) + |d S . 
Jo 

□ 

Remark 8. TTie SDE for the approximated is easier and faster to simulate 

than i3.1\) . as it is exhibited by the following example. Notice additionally that 
I/| J ' ei ' is a continuous functional of the process (|3.4[) and of the Brownian path 
Wf . Eventually, by using L[ l ' ei ^ as the LIBOR rates, the computational complexity 
of the drift and thus of the model can be reduced substantially, while maintaining 
accuracy of prices. 

Example 2. In this example, we examine the performance of the strong Taylor 
correction method. Let N = 3 and consider pricing a caplet on the LIBOR rate L 1 
with strike K . Its price is given by: 

P cpt = aE P ((L 1 Ti -K) 



Assume that the volatility functions o~ l {t) : [0, Tj] — > M for i = 1,2,3 are given by 
(cf. Brigo and Mercurio |BM01| , formulation (6.12)): 

a\t) = (a(Ti -t) + d) exp ( - 6(2} - t)) + e, 

where the constants a, b, d, e are the same for all three LIBOR rates and are equal 
to a = -0.113035, b = 0.22911, d = -a, e = 0.684784. Thus, we can write the 
model under the terminal measure P as: 

+ aX ( t 

■(2,ei) _ 2f,, r (2,e 1 ) f aXf^ a 3 {t)p 2 



V l + aXp ei > 1 + aXl' 1 ' ' 

dhf^ = a\t)L^ ( - aX t 3 ' £l)a y )p » 3 )dt + a\t)Lf^dWl 

\ 1 + aX (3, £l ) ) 



dL 3 = cr 3 (t)L 3 dW 3 , 

dX (hex) = ,J.i^y(Mi)/ aXl 2 ' ei) a 2 {t) Pl2 aX^' ei) a 3 (t)p 13 \^ t _i^ v (i, ei ). 



(o-\t)X^ ( - aAt " W f 2 - ,iT )dt + a\t)X^dWl 



dxf^ = ei 



V V 1+alf' 1 ' y y 



dX t 



(3,ei 



initial values L^'^ = Xq 1 ' 61 -* = c,, for i = 1,2,3 and for all e\ . The Brownian 
motion vector {W} ,Wf ,Wf) is correlated with correlation coefficient pij given by: 

Pij = 0.49+ (1 - 0.49) exp (-0.13|i - i,j = 1,2,3. 
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The SDEs for the approximated LIBOR rates t^' ei ^ and L[ 2 ' ei ^ are given by: 

dLt ~ {t)Lt V - l + a(c 2 + e 1 Y?) + l + a (c, + e l Y?) + r + 

+ <j l {t)L^ ) dWl 

dlf^ = 0*®%™ ( - ^f^^ + «\t)L?^dWl 
4 w * V l + a(c 3 + ei r f 3 ) + / w * * 

T/ie partial derivative terms and Y t 3 are equal to: 

Y t 2 = cJ [\^ S )W?-^- f t a 2 (s)aHs)ds), 

Wo 1 + ac 3 Jo ' 

Y t 3 =c 3 [ o- 3 (s)dW*. 
Jo 

We compare three caplet prices: 

• benchmark price, underlying L^ 1 ' 61 ^; 

• strong Taylor price, underlying l[ 1,£1 ^ ; 

• frozen drift price, underlying L^ ,0 \ 

Numerical results in basis points (bps) are displayed in Table [H for parameters 
ei = 1, N = 3, a = 0.50137, ci = 3.86777%, c 2 = 3.7574%, c 3 = 3.8631%, 
T\ = 1.53151, Ti = T\ +ia, i = 2,3,4. We characteristically observe the difference 
in prices between the benchmark and frozen drift price, whilst our strong Taylor 
correction method performs very well and is computationally simpler and faster. 



strikes 


K=3% 


K=3.5% 


K=4% 


K=5.75% 


K=6.25% 


K=8% 


benchmark 


11.1831 


8.5897 


6.5503 


3.0349 


2.4423 


1.2969 


strong Taylor 


11.0687 


8.5691 


6.5867 


3.1448 


2.5513 


1.3926 


frozen drift 


13.9551 


11.1822 


8.8803 


4.6313 


3.8506 


2.2524 



Table 1: Caplet values in bps for parameters ei = 1, a = 0.50137, c\ = 3.86777%, C2 = 3.7574%, 
c 3 = 3.8631% and Ti = 1.53151. 



3.3.2. Weak Taylor Approximation. In what follows, we provide some results on 
how to correct option prices obtained by the SDE with the frozen drift (|3.5|) by 
adding a correction term involving the appropriate Malliavin weight. Let L^ fc ' £l 

denote the vector of the LIBOR rates (Lj,' ei \ . . . ,L^ ei ^). 

Proposition 2. Assume the setup of Lemma[ll where the i th LIBOR rate is given 
by: 

N v-tf.ei) 

(3.8) dL? ei) = a\t)L^( - Y '-^± PlJ dt + dWl 



with Lq' £i ) = Xq' 61 ^ for all t\ and all i = 1, ...,N. Assume furthermore that the 
Malliavin covariance matrix 7(L^ fe ' ) is invertible. Then the price of an option with 
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payoff g{J2f' ei ), for i < k < N and g bounded measurable, can be approximated by 
the weak Taylor approximation of order one: 

(3.9) Wl(g,L^) = P(0,T)(E P ( 3 (L^°)) + e 1 E r (g(L%<°)( Ti 
where the Malliavin weight (j^ is given by: 



CUD) =.>((ALV / -*V: I 



£1=0 1 



/or* < 2*. 

Proof. The weight £7^ is obtained by (|2.4|) . Notice that we can write: 

E p ( 5 (L^ 1 ))=E p ( 5 (L^°)Ct ! ), 



dei 



£1=0 



and hence the result (|3.9p by Definition [5] for n = 1. □ 

Example 3. In i/iis example we let N = 3 and we price a payers swaption with 
strike price K and maturity T\ , where the underlying swap is entered at T\ and has 
payment dates T 2 and T 3 . We assume that the volatility functions o~ % (t) : [0, TJ — * K 
for i = 1,2,3 are constant: 

o- l (t) = ai,a 2 (t) = a 2 ,o- 3 (t) = 03, 

such that we obtain under the terminal measure P: 

, r (i,£i) r (i.£i)/ / aX t °~i aX y t V 3 \ x \ 
dL '=a x L\' '[Pi2{ 75— r w -r)dt + dWI), 

v-(3i«i) 

dLp' £l) - - r (2 ' £l) ^ - 1 ° z 



V 1 + aXp ei) ' 
(3.11) dL\ = a 3 L%dW?, 

ir(2,ei) ta(3.ei) 

^=e 1 (. 2 xf-)(--^^, t + ^)) ! 
=e 1 ( ( 7 3 l i (3 ' £l W ( 2 



uu'i/j initial values Lq'^ — Xq^ 1 ^ — ci, for i — 1,2,3 and for all t\. W} and 
Wf are correlated with correlation coefficient p\ 2 . We freeze the drifts in the above 
equations to obtain: 

r(i.o) ( wi ( 1 "C2CT2 . ac 3 a 3 . 1 \ 

L = Cl exp CTiiy t - (Pi2(— + — ) + -ai)a 1 t 

V V 1 + ac2 1 + ac 3 2 / 



r (2,0) / TJ/ 2 / ac 3 cr 3 1 \ 

i = c 2 exp [a 2 W t - I — + -CT 2 Iff 

V M + ac 3 2 / 

L? = c 3 exp (a 3 Wf - -aft) . 

Similarly to the previous example, we compare four option prices: 



• benchmark price; 

• frozen drift; 
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• strong Taylor price; 

• weak Taylor price. 

The weak correction formula (|3.9| adds a correction term to the closed form price 
of the option. The swaption payoff at Ti can be found for example in [MR98 : 

N+l fc-1 
k—i+1 j—i 

if the underlying swap is entered at time Ti and has payment dates Tj+i, T. oik 
is given by: 



Oik 



Ka, k = i + l,...,N, 
l + Ka, k = N+l. 



The payers swaption value at time t = can be written as: 

pswptn = p(0 7 T,)E P ,(P^ P * n ) = 
JV N 

(3.12) =P(0,T)E P ((-^a fc [](l + aL^)-(l + if a )) ), 

k—i j—k 

where cti := —1 and P 1 denotes the forward measure corresponding to the bond 
P(t, Ti) as numeraire. Therefore, its benchmark price is given by the above formula 
with N = 2 and i = 1 : 



/ + 



bp sw P tn = p ^ j,^ I jgp^Q,^ + aL ^ + qPl^l^ - Ka 2 L 2 Ti - 2Ka) 
Its weak Taylor price is given by (|3.9p with i = 1 and k = N = 2. 
wP° wptn = P(0,T) f%((aL£' 0) + aL ( *f + a 2 L^L$f - Ka 2 L^ - 2K 

+ ei E P ((a^; 0) + aL^ + c?L%?>I%f> ~ Ka 2 L$f - 2Ka) + ( Tl 

The weight £n is given by (|3 . 10[) . The partial derivative terms := ■^■\ ei= oL^ ei ' > 
and C\ x := -^-\ ei —nL^ ei ^ are given by: 

(i,o) f Tl _ ( V3ac 3 f3 2 

and. 



1 .In V l + acs) / 



o V(l + ac 3 ) 

2 _ r(2,0) /" Tl _a TJr/ 2, 



= L^> I ~o- 2 (3 3 W?dt, 



2 2 2 2 

with Cq — Cq — and /3 2 := i^^^i ' @ 3 : ~ (h+aca) 1 • ^ e Malliavin covariance 



matrix of the vector (L^°\ L^?' ^) is eguaZ £o: 



(l + p? 2 )(i&° ) ) 2 T 1 a? 2p 12 (4; 0) )(^ 0) )T 1 a 1( T 2 



7 ((Lg'°),Lg' ))) = 



2 Pl2 (i^ 0) )(i(?; 0) )T 1 a 1 a 2 (1 + ^(Z^) 2 !^ 
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The determinant is not zero as long as pu ^ 1, which is a natural assumption. 
Hence under this condition, its inverse is given by: 



7 -i((L£V£ 0) ))- — 



(1-P? 2 ) 



/ I+P12 2p 12 \ 



2pi2 I+P12 



Write the weight £ti = Cji + Cti > where the first weight is obtained as: 



and similarly: 



C 2 Tl = f (D 2 t 4f\c l Tll ^ + C 2 Tllu l ) + D 2 t L {2 f\c l Tll ^ 1 + C 2 Tll ^))6W? 



Performing all necessary calculations, we conclude that: 
Analogously we obtain £j. as: 



2 Ti 



Tl ' W 2 dt 



1 

Notice that the weights are functions of normal variables and thus the calculation of 
the weak Taylor price amounts just to computation of deterministic integrals. Table 
\3\ gives the swaption prices in bps for parameters N = 3, a — 0.25, o~\ — 18%, 
cr 2 = 15%, ct 3 = 12%, c = 5.28875%, a = 5.37375%, c 2 = 5.40%, c 3 = 5.40125% 
and p\2 — 0.75. 



strikes 


K=4% 


K=4.5% 


K=4.75% 


K=5% 


K=5.15% 


K= 


-5.25% 


benchmark 


10.2240 


6.5386 


4.7454 


3.1060 


2.2599 


1 


7758 


frozen drift 


10.2132 


6.5326 


4.7419 


3.1028 


2.2582 


1 


7618 


strong Taylor 


10.2240 


6.5386 


4.7454 


3.1060 


2.2599 


1 


7758 


weak Taylor 


10.2266 


6.54 07 


1 7485 


3.1064 


2.2593 


1 


7626 



Table 2: Swaption values in bps for parameters ei = 1, a = 0.25, a\ = 18%, <T2 = 15%, 
0-3 = 12%, c = 5.28875%, ci = 5.37375%, c 2 = 5.40%, c 3 = 5.40125% and m = 0.75. 
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3.4. The Stochastic Volatility LIBOR Market Model. In this section, we 
develop a stochastic volatility LMM. The stochastic volatility parameter vt follows 
a square root process, like in the extensively applied Heston model [Hes93 . The 
resulting model, called hereafter the stochastic volatility LMM (SVLMM), has the 
following dynamics under the terminal measure: 

(3.13) dL\ = cj i {t)L\^r t (- V a tCX ) Pij jv t dt + dWl),i = l,...,N, 

dv t = k(9 - v t )dt + e 2 \/vtdB u 

where k,8,€2 £ R+. The Brownian motions Wt = (W/, Wf*) and B t are ex- 
pressed under the terminal measure with correlations dW?dB t = pidt and dW^dW? = 
Pijdt for i, j = 1, ...N. We assume additionally that the filtration {Tt)o<t<T is gen- 
erated by both Brownian motions. Observe that the process v t is a time-changed 
squared Bessel process with dimension 8 = AnO/e^. If 5 > 2, then the point zero is 
unattainable. So we require 2k9 > e| for the process vt not to reach zero. 



3.4.1. Pricing a multi-LIBOR option. In this section, we aim at approximating the 
price of an option with payoff depending on the vector L^ fc ' ei:£2 = (L^ £l ' C2 , . . . ,L^' ei ' £2 ). 
We interpret the volatility of the volatility parameter ti as a parameter on which 
the LIBOR rates depend. Overall, we parametrise the SVLMM by both e\ and €2 
and correct prices in a weak sense introducing Malliavin weights. 

Proposition 3. Consider the SVLMM (|3 . 1 3[) and assume that the Malliavin co- 
variance matrix 7(L^,' ,0:0 ) is invertible. Then the price of an option with payoff 

7/>(L^' £l ' e2 ), i<k<N, where ip is a bounded measurable function, can be approx- 
imated by the weak Taylor approximation of order one: 

W^^L^ 2 )) = P(0,T)(E P (^(L^ o < )) +6 1 E p (^(L™)Ct 1 ) + 

(3.14) +e 2 E P (^(L™V Ti )), 

where the Malliavin weights £r< , ttt; are given by: 



(3.15) 
(3.16) 
fort < Ti. 



Ct, = 6{[D t Lj,. ) 7 (Ly,. ) — 

x(t r> T i,fc,0>0\T -1/T »,fc,0,CK 9 
7T Tj = ^(AL^ ) 7 ( L tV 



■ i.k.e 



ei =o 



rp 



Proof. The weights (i; and ^ are obtained by ()2.4|) . We derive (|3.14|) by noticing 
that: 



E(^(L^' £1 ' £2 )) 



= E(^(L^°))+e 1 A 



£1=0 



E^(L^'°))+e 2 A 
Ce2 



£2=0 



E(^(L^'°' £2 )) 



= E(V^L^ )) + £i E(^(L^°)CtJ + e2 E(^(L^°) 7 r Ti ), 
from Definition [2] for n = 1. 



□ 
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Example 4. Let N = 2 and consider the SVLMM where the volatility functions 
<r*(t) : [0, Ti] —y R for i = 1, 2 are assumed to be constant and in particular o~ l {t) — 
o~\, o 2 {t) = <J 2 . We derive an approximative formula for the price of a payers 
swaption with maturity T\ and strike price K . The underlying swap is entered at 
T\ and has payment dates T 2 , T3 . Under the terminal measure P we can write the 
SDEs for the LIBOR rates and stochastic volatility as: 

dv c t 2 =k(0- vf}dt + e 2 y/tf*dB t , 

, r (l,ei,e 2 ) r (l,ei,e 2 ) aA t CT 2 t2 ,, r (l,ei,e 2 ) ATJjjTrl 

dL\ = —L\ ' p\2 n ?(TlV t a dt + <TlLi 'Wv t 2 dW t , 

dL?^=* 2 Lf^J^dWl 



dX { t 2 ' ei) = e 1 {a 2 X^ t2) ^vt 2 dW i 

W{ and W 2 are assumed to be correlated, so correlations are as dW^dB t = p%dt 
and dW^dW 2 — p\ 2 for i = 1, 2. The (0, 0)-model is given by: 

u° =exp(-Kt)(v%-e)+0, 



r (1.0.0) ( f rojurl ( ac 2pl2 . 1 s 

LV = ciexp [ax I \ v^dW t - (— a 2 + -<Ji)cai 

l£ 0) = c 2 exp (a 2 ' ^fidW 2 ~\v\c 



v(2,0,0) _ 
A t - c 2 , 

c := f^ 1 i^dt = 0Ti — V ° K (exp (— kT±) — 1). ^4s m the previous example, we 
compare the following option prices: 

• benchmark price; 

• frozen drift; 

• weak Taylor price f|3 . 14[) . 

The benchmark price is given by (I3.12p wii/i N = 2 and i = 1: 

bP ™Ptn = p^ T ) Ep ^ aL (^Ut2) +aL ^ 
The weak Taylor price is obtained by (|3.14[) ; 

w ps Wp tn = ^^ aL U,0,0) + a£ g,0) + a 2 L {lfi,0) L ^fi) _ Ka 2 L (2fi)_ 

- 2^a) + ) + ei E P ((aZ£ A0) + aig' 0) + a 2 L^°' 0) L^ 0) - i^a 2 - 
• ig'°> - 2Jfa) + C Tl ) + e 2 E P ((aL% m + aL™ + a 2 L%°> 0) L$f- 



Ka 2 L {2 f ] - 2Ka) + TT Ti y\. 



We calculate the Malliavin weights Qt-l , ttti o.s given by (|3 . 1 5[) and (|3.16[) corre- 
spondingly. We can express the weight Qt-l o,s: 

Cti = Ct\ + Cti > 
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with: 



d 
de 1 



£i=0 



L (l>ei,0) 



-l/ r (l,0,0) .(2,0,0), \ 
[L Ti ,L Ti juj + 



/ n w (2,0,0) ^ 



£1=0 



- i (4 i ^,^) 21 )U^, 



and: 



r (l,ei,0) 1^(1,0,0) r (2,0,0)N \ 



( Tl -J o \[D t L Ti — 



swf. 



since ^-| ei=0 L^ €l ' ^ = 0. The partial derivative term with respect to e\ for L 1 
given by: 



_d_ 



£1=0 



r (l,€l,0) _ j(l 



rTi 

L,0,0) / o/ 

1 7o 



v°aW 2 )di 



(1 + «C 2 )2 7 Q 

= -o- lPl2 p2L% m J 1 (fl(Ti - t) - ^^(cxp (-kTi) - exp {-Kt)))dW 2 
where f3 2 — ci+acty 1 ' Similarly the weight tt^ is given by: 

1 i 2 
7T Tl = 7T Ti + 7T Ti , 



with: 



j ° i=i j=i 



d_ 

de 2 



£2=0 



^ 0,ea) 7- 1 ((^ 0,0) ,^° ,0) )«))*Wi 1 . 



and: 



Partial derivative terms are equal to: 



_d_ 

de 2 



r (i,o,£ 2 ) r (i,o,o)/ffi [ 1 exp(-Kt) f rn, n „ i 

=° v z Jo Wvy Jo 



+ 



1 a-f t ac 2 a 1 a 2 p 12 . 



) I exp (ks) -\/^° ( exp (— k7\) — cxp (— ks 



k 2 1 + ac2 
Doing similar calculations, we derive the second partial derivative: 



d_ 

de 2 



£ 2 =o 



r (2,e 2 ) _ r (2,0)/l 



v u t " / ° 



w/iere Vt = cxp (— «t) J* exp (Ks)\/v®dB s . 
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We calculate the Malliavin covariance matrix 7((Ly i '°'°\ L^ 0,0 -*)) and its in- 



verse. 



(, 



(l + p? 2 )(L^°' 0) )V / v»dt 2p 12 L^"L^ , 'a 1 a 2 c 



(1,0,0) j- (2,0,0) 



\ 



V 2p 12 4Y'° ) L^a 1 a 2 c (1 + pl 2 )(L^)^ 2 c J 

^det( 1 ((4f-°\^ 
Hence its inverse is given by, for p\ 2 ^ 1: 



1 



(1-Pf 2 ) 



V 



1+PI2 

(^•°' 0, )V?c 



2pi2 
r (1,0,0) , (2,0,0) ~ 



2pi2 



, (1,0,0) r (2,0,0) 



(L<? ; '°))V|c / 

If we define Xi = J^ 1 y/v°dW t \ i = 1,2 and Y = J^ 1 ^(^-^-^(exp (-kTi)- 
cxp (— nt))jdWt, we finally obtain the weights as: 

and: 

C 2 Tl =^(x 2 Y-Cov(X 2 ,Y)). 
Moreover, for the weight tt^ we define: 

r Ti 



B = L ff W( cxp ^ KTl ) -exp(-/ct))dB t , 



and random variables Di, Zi for i = 1,2: 



A 



Tl /(*) f g{s)dWidWl 

Jo 

Ti 



^ = / 1 /(*) / sWdZjdW?, 
Jo Jo 



where the Brownian motions Z\ are independent from W{ and f(t) = 22lLjS*2_ 



g(s) — exp (ks) . Therefore, we obtain the weights as: 

1 1 / , /: T v B^i ^c 2 (2p 12 g 2 + 
1 ~ 2c [ X ^ D ^ + V 1 - + — \ 1 + J " 

P\E ( ac 2 (2p 12 o 2 +cti) + p i2 



K V 1 + «C2 







2c 



Xi(p 2 ^ 2 + A/l-p^ 2 ) + 



+ — — 1 + 

K AC 
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where E equals to ( ~ (l — exp (— kTi)) — T\ 

( 



1 

2c 



(7 2 



X 2 {p2D 2 



P 2 2 Z 2 ) 



D 



Similarly we get tt^ as: 
, <r 2 p 2 E\ 



2c \ cr 2 



+ 



P12.B 



1 - p\Zi) - 
aip 2 E fac 2 (2pi 2 a 2 + a 1 



-{o 2 X 2 + ] 

o\BX 2 f ac 2 (2pi 2 a 2 + a\) + o\ 
1 + ac 2 



na 2 



1 + ac 2 




In this example, the weights are functions of normal variables and double sto- 
chastic integrals, which are computed via simulation. Table ^ reports the swaption 
prices in bps with parameters N = 2, a = 1.5, <7i = 25%, a 2 — 15%, cq — 5.28875%, 
ci = 5.4%, c 2 = 5.39%, v = 1, p! = -0.75, p 2 = -0.6, k = 2.3767, = 0.2143, 
e 2 = 25%, pi 2 = 0.63. 



strikes 


K=3.5% 


K=4% 


K=5% 


K=6% 


K=7% 


K=8% 


benchmark 


3.8984 


2.9221 


1.2588 


0.3858 


0.1019 


0.0216 


(0, 0)-model 


3.8951 


2.9053 


1.2705 


0.3966 


0.0942 


0.0185 


weak Taylor 


3.8990 


2.9159 


1.2694 


0.3791 


0.1042 


0.0210 



Table 3: Stochastic volatility swaption values in bps for parameters ei = 1, a = 1.5, a\ = 25%, 
0-2 = 15%, c = 5.28875%, ci = 5.4%, c 2 = 5.39%, -uo = 1, pi = -0.75, p 2 = -0.6, 
ft = 2.3767, 9 = 0.2143, e 2 = 25%, pi 2 = 0.63. 
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